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Staggered moment dependence on field-tuned quantum fluctuations in 2D frustrated 

antiferromagnets 
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We propose an efficient method to identify the degree of frustration in quasi-2D antiferromagnets described 
by the Ji — J2 Heisenberg model. The frustration ratio J2 /J\ is usually obtained from analysis of susceptibility, 
specific heat and saturation field. We show that the non-monotonic field dependence of the staggered moment 
caused by the suppression of quantum fluctuations in a field depends strongly on the frustration ratio. This gives 
a powerful criterion to determine J2/J1 using a combination of exact diagonalization (ED) method for finite 
clusters and spin wave analysis. We apply this method to the quasi-2D compound Cu(pz)2 (0104)2 and show 
that it leads to an intermediate ratio J2/J1 — 0.2 for the frustration. We also explain the observed anomalous 
increase of transition temperature in applied fields as an effect of reduced quantum fluctuations. 

PACS numbers: 75.10.Jm, 75.30.Cr, 75.30.Ds 



I. INTRODUCTION 

Quasi-two dimensional antiferromagnets may show a sub- 
tle interplay between quantum fluctuations and exchange frus- 
tration. This is clearly apparent in the 2D square lattice Ji — J2 
model which has been found to be approximately realized in 
a class of layered (S — 5) compounds'-^. An anisotropic 
version of this model is also relevant for Fe pnictides^^. The 
ground state of the model is determined by the size of the frus- 
tration angle (j) ~ tan^^ {J2/J1) (—7Z<(j)<7t). Depending on 
(j), Neel as well as columnar antiferromagnetic (NAF/CAF) 
structures with wave vectors Q={n,n) or (7r,0) respectively 
may be realized. The size of staggered moment in these 
phases is reduced from the classical value 1/2 by the effect of 
quantum fluctuations in the ground state. It depends strongly 
on the frustration angle and vanishes around ^/tt ~ 0.15 
(J2/J1 = 0.5) where a nonmagnetic stacked dimer state'* and 
close to 0/;r ~ 0.85 (J2/J1 = —0.5) where a spin nematic 
ground state^ appears. Therefore the determination of is 
of central importance to characterize frustrated 2D quantum 
magnets. Frequently it is obtained by comparing results of 
high temperature series expansion or finite temperature Lanc- 
zos method (FTLM) to the temperature dependence of sus- 
ceptibility and magnetic specific heat^-^. Analysis of the sat- 
uration fields is also employed^. However, in contrast to the 
determination of a single J for the nonfrustrated model, the 
former method is inaccurate and even ambiguous^-^ for the 
frustrated case. It is not able to discriminate between frustra- 
tion angles and j — (p, in particular this means if J2/Ji< 
the cases Ji > 0,^2 < and 7i < 0,^2 > cannot be distin- 
guished. 

Here we propose a powerful method for the determina- 
tion of the frustration ratio in 2D quantum antiferromagnets 
in each sector. Quantum fluctuations depend on the cant- 
ing angle between the moments, it was shown in Ref.^ that 
this leads to a nonlinear uniform magnetization characteris- 
tic for the degree of frustration. An even more drastic effect 
may occur in the staggered moment'". The classical canting 
together with field suppressed quantum fluctuation lead to a 
non-monotonic field dependence which depends crucially on 
(j). Comparison of theoretical prediction from numerical ED 



and spin wave theory with experimental results of the stag- 
gered moment field dependence mQ{H) can precisely deter- 
mine the frustration ratio in the NAF and CAF sector. As an 
example this method is demonstrated for the 2D Cu-pyrazine 
compound Cu(pz)2(C104)2. We show that a fit to mQ{H) and 
the FTLM fit to x{T) both give values for the frustration angle 
(p which are considerably larger than reported previously al- 
though they are still deep within the NAF region. We demon- 
strate that the former method gives a reliable value for (j) in 
contrast to the latter. In addition the anomalous field depen- 
dence of the transition temperature Tn{H) is explained within 
a self-consistent RPA spin wave theory as a result of frustra- 
tion suppressed by the field. 



II. STRUCTURE FACTOR OF THE J2 MODEL 

Our analysis is based on the quasi-2D frustrated square lat- 
tice spin S=l/2 model including the Zeeman term: 

— JlY^Si -Sj +J2 52 ^i'^j 

(U) (('■;■» 
+f Y^SrSj-g^sHY^S, (1) 

Here Ji — Jc cos and J2 = Jc sin are the in-plane near- 
est and next nearest neighbor exchange constants with Jc = 

{jI + J2)- giving the overall energy scale. Furthermore J' 
is the coupling between 2D layers with J' /Jc <C L The lat- 
ter term will only be needed for analysis of 3D AF transition 
temperature. The field dependent total moment is given by 



mJ{H)=ml{H)+m}^{H) 



(2) 



consisting of uniform and staggered moment perpendicular 
and parallel to the plane respectively (inset of Fig. [Hi. Within 
the numerical ED Lanczos approach they can be expressed 
in terms of the static structure factors by = Svv{Q) and 



\Sxx{Q) + Syy{Q)\ using the definition 



1 ^ 



(3) 
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FIG. 1 . Field dependence of staggered moment from spin wave the- 
ory for several frustration angles <j>. When <j) approaches the (classi- 
cal) critical value <p /k — OAS where AFM at zero field breaks down, 
»jq(/i) behaves strongly non-monotonic. 



The expectation values should be evaluated with the corre- 
sponding ground-state at field H. The uniform and staggered 
moments are obtained by scaling the structure factor with the 
size of the clusters. At // = 0, the spin (S) dependent nor- 
malization factor is -yVap =N{N+ l/S)as explained in Ref.|3|. 
It must be modified for finite magnetic fields to ,y^^j. = A^^ and 
.yVicx = = N{N — 1) to account for the effect of dropping 
non-zero on-site terms for the staggered moment in the above 
sum, Eq. OJ, which is necessary to achieve the correct limit- 
ing values at the saturation field for arbitrary tile size A^. 



III. UNIFORM AND STAGGERED MOMENT FROM 
LINEAR SPIN WAVE THEORY AND EXACT 
DIAGONALIZATION 

It is instructive to calculate both moments also within spin 
wave approximation for comparison. In the classical limit 
one has (in units of gjiB)' '«o ~ Scos&c, '«Q ~ Ssin&c 
where cos 0c- = H/Hs is the classical canting angle of spins 
as counted from the field oriented || a or perpendicular to the 
quasi-2D be plane . We use this convention to be compati- 
ble with Cu(pz)2(C104)2 discussed later For the NAF case 
the saturation field is given by Hs — 4Ji/{gHB) which is in- 
dependent of J2- In the classical picture the staggered mo- 
ments of constant size (H) — S aie simply tilted out of the 
plane until at only the uniform moment niQ = S is left (in- 
set of Fig. [Til. For H = quantum fluctuations reduce the size 
of the moment and the amount of reduction strongly depends 
on the frustration degree given by 0. However, at saturation 
(H — Hs) the total moment always has to return to the clas- 
sical value nif ^ niQ ^ S because the fully polarized (ferro- 
magnetic) state is an eigenstate and quantum fluctuations are 
absent. This means that the field dependence of wq(/!), mo(/i) 
and likewise of mt{H) is tuned by 0. Here mo{h) is the uni- 
form moment or magnetization which was shown in Ref.'o' to 
exhibit nonlinear field variation depending on the size of (j). 
Similarly the field variation of the staggered moment mQ{h) 
should be strongly influenced by the size of frustration. We 




FIG. 2. Parametric representation of uniform (ino{h)) and staggered 
(mQ{h)) moment as function of field. Here /f = at the lower right 
and H = Hs at the upper left comer. Dashed line: classical result. 
Symbols: ED values for finite tiles with = 8 — 30 and (j)/7t = 0.063. 
Empty (full) circles: extrapolated values from finite size scaling for 
(l>/7t = 0.04 (<p/7t = 0.063). Full thin (thick) fine: spin wave result 
for intermediate frustration (jy/n = 0.04 (0/;r = 0.063). Dash-dotted 
line: spin wave result for small frustration (j>/7t = 0.006. 



propose that this effect may be used as powerful means to de- 
termine = tan^^ (72/^1)- 

The Hamiltonian ([T]) may be diagonalized in the local coor- 
dinate system of the canted spins^ leading to spin wave ener- 
gies 



Bk=Jk--iJ} 



BkCos2 0,]2-[Bk(l 
A-q) 



cos^ 0, 



1 , 

2 (A 
1 



(4) 



k+Q 



where 7k ~ 2Ji {cos ky + cos k^) + 4^2 cos A:,, cos fc^ is the 
Fourier transform of the intra- layer exchange couplings and 
k and Q are wave vectors lying in the bc-plane. Here we 
set J' — 0. Because f /J ^ 1 it is not important for ground 
state properties. Defining h = gpLsH we have cos 0c = h/hs 
with hs — 2SAq. The uniform and staggered moments are then 
given by the ground-state expectation value of the total mo- 
ment projected onto the field direction and the plane perpen- 
dicular to the field, respectively, leading to 



mo — cos 0c 



^ 1 1 y. gk(Ak-gk) 



(5) 



sin 0c < 5 



1 



2sin^ 0c 



1 ^Ak— Bkcos^0c 



N 



Ek{h) 



cos2 0c l y.Bk(Ak-Bk) 



Isin^&rN' 



(6) 



up to ff{l/S). The calculated staggered moment for various 
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FIG. 3. Uniform magnetic susceptibility from FTLM for various tile 
sizes. Fitting to experimental values (dots) for Cu(pz)2(C104)2 starts 
above the vertical line and below the maximum. Optimal fit parame- 
ters J\ , J2 are indicated. The g- value is obtained from gjlsHs = 47i . 



frustration angles is shown in Fig. [T] It is obvious from 
this figure that the field dependence of the staggered moment 
mQ{h) is strongly influenced by the frustration angle. In con- 
trast to the classical case it exhibits non-monotonic behavior 
This appears because firstly the size of the total moment in- 
creases with field due to suppression of quantum fluctuation 
and secondly the moments are canted out of the be plane, 
which reduces the staggered projection. For larger (j) one ap- 
proaches the region of the nonmagnetic phase ((p/n ~ 0.15) 
where spin wave theory eventually breaks down (mQ(0) — !• 0). 
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FIG. 4. Field dependence of staggered moment squared (~ scat- 
tering intensity) in Cu(pz)2(C104)2(full circles) (from Ref."). Fat, 
thin and dotted lines: Results from linear spin wave theory (LSW) for 
(I)/tz = 0.063,0.04,0.006 respectively (y2/Jl= 0.2, 0.126, 0.02). Dia- 
monds: ED results for (^/tt = 0.063. Top: Scaled up to the available 
experimental values. Bottom: Scaled up to the saturation field. 



We have checked the results of spin wave theory with an 
unbiased numerical ED approach for finite size clusters. The 
uniform and staggered moment are obtained from the scaling 
analysis of the calculated structure factor performed in anal- 
ogy to the zero-field case^. The results are shown in Fig. |2] 
in a parametric representation where m(j{h) is plotted versus 
mQ{h). The field is the parameter in the plot starting at H=0 
in the lower right corner (zero field staggered moment) up to 
H = Hs in the upper left corner (saturated uniform moment). 
The classical reference curve is shown by the dashed line. 
The (unsealed) moment values for various sized 2D J\ — J2 
tiles with intermediate frustration angle ^/n = 0.063 are pre- 
sented as symbols. The finite size scaling extrapolation to the 
thermodynamic limit has to be performed for constant mag- 
netization given by the horizontal dashed lines. The result 
is shown in open and full circles. They agree with the pre- 
dictions of spin wave calculations (full lines) for (p/n = 0.04 
and 0.063. For comparison the spin wave result for the the 
nonfrustrated (0 = 0) NAF (dotted line) is also shown. It is 
suggestive from Figs. [T]and|2]that a careful determination of 
the ordered moment field dependence may determine the de- 
gree of frustration in a quasi-2D antiferromagnet given by (j) 
ovJ2/Ji. 



IV. APPLICATION TO QUASI-2D Cu(pz)2(C104)2 

To demonstrate the strength of this method we apply it to 
the quasi-2D antiferromagnet Cu(pz)2(C104)2- In previous 
work'^'^^ the spin wave excitations, field dependent moment 
and ordering temperature Tn{H) of this compound were deter- 
mined. Based on these results and on earlier thermodynamic 
studies'^''^ it was proposed that Cu pyrazine is close to a pure 
NAF with a frustration ratio ^2/^1 = 0.02 ov k 0.006. 
The field dependence of staggered moment however was not 
analyzed in this respect. This will be performed in the present 
work to illustrate this powerful method for obtaining Ji/Jx- 
We also compare to the results of the analysis of x{T) us- 
ing our FTLM data. The previous work^^'^^ based on series 
expansion assumed from the outset that 72=0 (0 = 0). We 
perform an unbiased analysis with possibly non-zero J2 by fit- 
ting the FTLM data with variable and Jc to the experimental 
data. In order to reduce the influence of finite size effects only 
data points from slightly below the maximum in x{T) up to 
the highest temperature are included. The result for the best 
FTLM fit from various cluster sizes is shown in Fig. |2] and 
Ji, J2 are given in the caption. They lead to a best fit with 
J2/J1 = 0.12 or (j)/7Z — 0.04. From the value of Ji and the 
measured saturation field Hj = 52 T we get a gyromagnetic 
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FIG. 5. Field dependence of transition temperature from RPA spin 
wave theory. Symbols: experimental values " . Full line: Calculated 
fromEq.l|7llwith0/;r = O.O63,//yi =6-10"^ Dashed line: (p/Tt^ 
0.006, J'/Ji =6.8-10-4 



ratio g = 2.11. For (p/K > 0.04, the agreement around the 
maximum of ;i^(r) becomes worse. 

However, the fit to x{T) does not necessarily yield an accu- 
rate and unique solution, since it depends on ©cw = -^i + Ji, 
but only weakly on the individual exchange constants and 
the complementary values (/) and 0' = — having identical 
©cw corresponding to Neel and columnar AF order cannot 
be distinguished from an analysis of x{T)- It is therefore im- 
portant to check this value of by using the field dependence 
of the ordered moment. The square of the latter is propor- 
tional to the scattering intensity. We calculated mQ{h) from 
the spin wave theory in Eqs. (|4]i and ^ using three different 
values of (j)/7Z. The results are shown in Fig. |4] The stag- 
gered moment squared (wjq ^ intensity) increases by about a 
factor of two in the measured regime up to H/Hs ~ 0.25. Us- 
ing spin wave calculation (which is accurate for ^/n < 0.1) 
supplemented by ED (see also Fig. |2]l the value (p/n ^ 0.063 
{Ji/Ji = 0.2) (full line) gives perfect agreement with experi- 
mental moments. This is somewhat larger than the value from 
the FTLM fit to xi^)- The experimental data in Fig. |4]are 
determined with 7 — 8% relative accuracyiJ-. 

For comparison we also show the moment for (p/n — 0.04 
(thin line) and the nearly nonfrustrated NAF with J2/J1 — 0.02 
01 (p/n = 0.006 (dotted line). In the latter model the predicted 
field induced staggered moment increase is much too small. 
To get a more pronounced moment increase with field one has 
to increase the frustration 0, as is evident from Fig. [T] and 
for (p/n = 0.063 the experimentally observed increase is ob- 
tained in Fig.m The discrepancy to the FTLM value of (j) may 
possibly be caused by the background subtraction process^' 
which gives an uncertainty to the the absolute size of moment 
increase. 

Quantum fluctuations not only lead to the distinct field de- 
pendence of the ground state staggered moment but are also 
responsible for the observed anomalous increase of the Neel 
temperature with field strength The Neel order at finite 7]v 
is an effect of the finite interlayer coupling J'. Because of the 
quasi-long range order of 2D HAF with exponentially increas- 
ing correlation length a small interlayer coupling J' /J ^ 1 
leads to sizable Tn on the scale of the intra-layer exchange 



strength'^. 7)v(0) for the nonfrustrated model (0 = 0) may 
be obtained from an empirical formula based on MC simu- 
lations, however, this is not available for <p ^0 and for finite 
fields. Therefore we use a self-consistent RPA theory based on 
quasi-2D spin waves where T]\i{h) is determined by the condi- 
tion of vanishing staggered moment. We obtain 



TN{h) 



N' 



Ak-BkCos-0,. 



{h) 



H -1 



(7) 



Now the interlayer coupling along a-direction is included 
according to 7k — 2J\{cosky + cos ^7) + 4^2 cos ^3, cos + 
2J' cosk^ and k as well as the ordering vector Q in Ak,Bk of 
Eq. (|4]l have now all three components. The integral in Eq. (|7]l 
is therefore finite leading to a nonzero T^ih)- Eq. (|7]i reduces 
to the expression in Ref . 16^ for h=Q. The RPA theory predicts 
the right dependence of 7)v(0) on J' /Ji but the absolute values 
are larger than from those of empirical formulas fitting the MC 
simulations'^. Therefore in Fig.|5]we have plotted the normal- 
ized field dependent transition temperature TM{h) /7iv(0). It is 
shown for two sets of values i(j),J'). The increase in 7iv(/i) is 
driven by the reduction of quantum fluctuations since in a field 
the average spin wave energy is increasing. The (p,J' values 
for the dashed curve reproduce the experimental Ta?(0) — 4.2 
K but fail for the field dependence. The values derived here 
(corresponding to full line) lead to excellent agreement with 
experimental Tff(hyTN(0) but 7/^(0) is about twice the experi- 
mental value. We think that the functional dependence Tn{H) 
is more significant for J' /Ji than the single value of Tn{Q). 



V. CONCLUSION 

We have presented the analysis of field dependence of or- 
dered moment in frustrated quasi-2D antiferromagnets using 
exact diagonalization for finite clusters and compared with 
spin wave theory. The staggered moment exhibits pronounced 
nonmonotonic behavior as function of field which depends on 
the degree of frustration given hy (p = tan^' [Ji/Jx)- This pro- 
vides a powerful means to extract the frustration ratio which 
is more accurate and less ambiguous than using temperature 
dependence thermodynamic quantities. We have applied this 
method to Cu(pz)2(C104)2 and conclude, primarily from the 
ordered moment field dependence, that it is a quasi-2D antifer- 
romagnet with intermediate frustration. This method may be 
used more generally for frustrated antifeiTomagnets. In par- 
ticular it should also be applicable when the field dependent 
ordered moment is extracted from analysis of NMR splittings 
rather than from neutron diffraction data. Since it requires the 
existence of an ordered moment it will, however, not be useful 
for compounds corresponding to the disordered regimes of the 
J\ — J2 phase diagrams, if they should indeed exist. 
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